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ABSTRACT 

A formulation is given for a collection of phonons (sound) in a fluid at a non-zero 
temperature which uses the simple harmonic oscillator twice; one to give a stochastic 
thermal “noise” process and the other which generates a coherent Glauber state of phonons. 
Simple thermodynamic observables are calculated and the acoustic two point function, 
“contrast” is presented. The role of “coherence” in an equilibrium system is clarified by 
these results and the simple harmonic oscillator is a key structure in both the formulation 
and the calculations. 
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1. Introduction 


The problem of understanding the thermal properties of a radiation field in a finite 
volume is both old and subtle. 1-3 Here a sound wave propagating in a water will be studied 
and the key issue will be the interaction of the sound radiation with the fluid matter and 
with the walls of the container. The time scales of sound waves, v — 20 — 2 x 10 9 Hz, 
and those of the water molecules lead to adiabatic (isentropic, if approximately reversible) 
thermodynamic processes rather than constant temperature transitions. 4 This work is a 
special case of a project by one of us (AVN) which addresses the full nonlinear problem of 
bubble formation by sound waves. The linearized problem will be studied here, where the 
fluid has a coherent interaction with the sound radiation and an incoherent, or stochastic, 
interaction with the reservoir. 

The harmonic oscillator has played a central role in the coherent states 6 11 and will 
be used for the coherent (Poisson) process describing the phonon radiation the sound field. 
Since the reservoir is incoherent the total interaction with the fluid is partially coherent. 14,15 
The reservoir is analogous to Feynman’s rest of the universe 12 and Han, Kim and Noz 13 
have shown the relation of this idea to quantum squeezed states and time-uncertainty. 

The model presented here will use the simple-harmonic-oscillator twice: first to gener- 
ate “stochastic or chaotic” noise and second to generate a Glauber coherent state of scalar, 
longitudinal phonons. This is a more realistic model of noise, in that it has both coherent 
and a random components. It will be called partially coherent following a standard useage 
in quantum optics. The density, entropy free energy and a two-point function which gives 
the acoustic contrast are all calculated. The reason that the SHO is so useful is that since 
their Gaussian functions are dense in L 2 , quantum mechanics guarantees that the crucial 
interaction between the fluid and the sound radiation can be approximated by an infinite 
collection of oscillators. This is why the approach of Planck 2 was correct even though 
quantum mechanics was not yet created. Also, finite energy classical solutions will lie in 
L 2 or at least in the Soboler space H 1 = L 2 ’ 1 , which is the space where the “function” 
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and its “gradient” are square integrable. 

In Sec. 2 the model will be presented, and the density p, the entropy 5, the free 
energy F and the pair correlation function are calculated for both single and N-mode 
partially coherent states. In Sec. 3 the Conclusions and Outlook are presented. 

2. The Model 

In Fig. 1, a schematic is given which shows a source of sound So (treated as a coherent 
state of phonons) a fluid F in thermal contact with the reservoir R, which is much larger 
than So or F. In general, phonons can enter the fluid from S'o and the fluid and reservoir 
can exchange particles as well as heat but all other exchanges are negligible. 



Fig. 1 . Schematic of the system modeled. The source of a coherent state of phonons 
is S'o, F is the stationary fluid volume and R is the reservoir which is much larger than the 
sum of S'o and F. The wavy lines indicate boundaries which allow particles and energy to 
pass. 
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The idea is a modification of one due to Kaup 16 ’ 17 that cavitons (here bubbles) are 
solitons (here solitary waves). In other cases, Williamsson and Wieland 18 , Glimm 19 , and 
others have shown that many physically interesting model solutions for plasmas and clas- 
sical fluids are nonlinear, coherent exitations of the medium. The formulation given can 
easily be generalized to Mi independent random components of noise and M 2 indepen- 
dent, coherent components. Thus, solitons and some other nonlinear modes could easily 
be added to the analysis. 

The phonons are bosons so that their creation and destruction operators a*, a satisfy 
the canonical commutation relations, 

[a, a] = 0 = [a*, a*], [a, a*] = 1 (1) 

and a unique, translationally invariant Fock vacuum |0 > exist s.t. 

a( x , t)\0 >= 0 . (2) 

The astersik power of an operator is its adjoint (a* and a are not self-adjoint) and on a 
complex number is its complex conjugate. Physically, the Fock vacuum is a quantum state 
with no phonons. The number operator N is defined as 


N = a* a 


( 3 ) 


and a number or Fock states is given by 


n >= 


(<■*)" 

77 .! 


10 > 


for each neZ + , the positive integers including zero. They are eigenfunctions of the number 
operator with eigenvalues nez+. The Fock representation Hf of the quantum Hilbert space 
'H is the L 2 closure of the linear span of the |n > states. The inner-product of the Hilbert 
space will be written as < -,- > and the inner-product compatible norm is written as 
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II . II = [< >]s. For any complex valued zeC 1 the unitary displacement operator, U(z), 

acting on the Fock vacuum yields the minimum uncertainty coherent state | z > given by 

oo n 

|,>=t?(,)|0>=e-W ! / j y:4=|n> . (4) 

n=0 v 

The Fock vacuum is the ground state of the SHO for the minimum uncertainty coherent 
states which will be used here. In terms of c-number coordinate q and momentum p the 
complex number 2 is written as 

z = (q,p) = q + ip ( 5 ) 


so that C 1 corresponds to the phase-space of the (1 — d) system. Two properties of the 
coherent states are that 


and 


a|z >= z\z > 


< x u z2 >= e - s ( * r_ ^ )(zi_i 2 ) e-* ( ^ Z2_2:iZ * ) 


( 6 ) 

( 7 ) 


From eq. (7) it is clear that the coherent are continuous in the label z and therefore are 
an overcomplete family of states, OFS. The L 2 — closure of the linear span of the coherent 
states |z > provides a continuous representation of the physical Hilbert space which will 
be written as 7f C s- 

A density operator is a positive, self-adjoint operator which satisfies 


P 


2 



( 8 ) 


The expected value of an observable A = A* in a state ipeTi with corresponding density 
operator p ^ can be expressed as 

< A >rp=< Vb Atf? >= Tr(p t pA) . (9) 

And additive thermal noise can be added “by hand.” The entropy of the system, 5*, is 
given by 

S = -k B Tr[plnp] (10) 
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where kg is the Boltzman constant. In information theory, one can set kg = l/ln2 and still 
use eq. (10). The entropy is obtained from maximizing eq. (10) subject to the constraints 


Tr(p ) = 1 


( 11 ) 


and 


Tr(pN) = c 


( 12 ) 


where ceR 1 is a parameter which labels the strength of the thermal state. If the thermal 
noise is Gaussian, its density operator can be expressed as 


p(c) = — f d 2 ze-^ /c \z >< z | 
nc J 


(13) 


In Rf this can be re-written as 


i r 00 2n 

p(c) = — j d 2 ze~\ z \ ! c e~\ z \ • >< ' n l 

7rc ^ n=0 


(14) 


Using Fubini’s theorem to interchange the integral and the infinite sum and then expressing 
z in plane polar coordinates gives 


M = £ 


•| n >< n| 


(c+lXl + l/c)^ 1 
in the Fock representation. A similar calculation of the entropy gives 


(15) 


S = kg[(c + l)/n(c + 1) — c /n(c)] 


(16) 


This sort of calculation was given by Glauber 10 , Wolf 14 , Sudarshan 15 and probably by 
others. The n t,l -order correlation function, gd n )(_XF ■ • • , -A n ), with Glauber’s normalization 
convention is 

y V \ 

9 (Xi, . . . ,X n ) 2n 


n [ca.»(x*,x»)] 

Jt=l 


1/2 
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(17) 


where the G^'^'s are Green’s functions or correlation functions. This is a quantum gen- 
eralization of the classical coherence degree 7 of Born and Wolf. 20 The visibility v of a 
two-slit interference pattern is related to 712 by 


v = l± |7i2(p)| 


(18) 


which physically represents the extremes in intensity. For the two-point function, </ 2 ^( - ), 
the G^ nm )’s are chosen s.t 


m Tr[pN(N — 1)] 

G {p) - rr(„AT)] J ' 


(19) 


This object is proportional to Glauber’s g ^ in ref. (10) but is not equal to his function 
because of different normalizations.. 

In the case of the thermal states, the Fock representation of g ^ for n ^ 0 is given by 


9 f \ l n >< n l ) = 


Tr 


\n >< n\N(N - 1) 




Tr(|n >< n\N) 


n * 


n 


( 20 ) 


A coherent state is very different from eq. (20) since 


Tr 


9 { S (I* >< A ) = 


\z ><z\N(N-l) 


Tr(|n >< n\N) 


- 1=0 , 


( 21 ) 


i.e. the Glauber state is perfectly coherent. 

For the thermal states, not surprisingly, g ^ has the opposite behavior from eq. (21) 
because 

1 \ n + ! | 


Tr 




E(*)(T4^) l">< nN(N-l) 


{rr[E(^T)(T^)” +, |n>< n |Jv]} 


E[(^l)(T+T7;)” +ln ( n_ ')] 


E(774TT)(T+r7c) n+1 ” 


T 2 


1 . 


( 22 ) 


(c+i) A 1+1/ 
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Defining 


eq. (22) becomes 


fn(c) 


(c + l)(l + l/c) 


n+1 


- JL 


E fn(c)(n 2 - n) - E fn(c)n 


tiL(c) = 


T 2 


= 1 , 


(23) 


E/n( c )« 

. n 

which is perfectly incoherent Now a partially coherent state has a density operator which 
is given by 


p(z;c) = U(z) PTh {c)U\z) = j- c Jd- 


■2 xe -\*\*/c 


\z + X >< Z + X | 


(24) 


A calculation similar to the one outlined between eqs. (14) and (15) gives 
Tr[ p (z-,c)N)=Tr{±- i <P X '-WK >< „|JV/ mn ] 

= |z| 2 + c . (25) 

It is straightfoward to calculate the entropy of the partially coherent state from 


S(z; c ) — —kgTr 


p(z;c)ln(pz-,c )) 


p{c)ln(p( c )) 


— —kgTr | 

= S(c) 

= kg [(c + l)/n(c + 1) - c /n(c)] , (26) 

where eq, (16) was used in the last two equalities. The pair correlation function, g i2 ' 1 , of a 
partially coherent state is given by 

12 \ 2 


9 ?} [*>(*; <0] - 1 ( c |j z p) ■ 


c + \z\ 2 .' ’ (27) 

in the coherent state basis. Let m > n with m,n€Z+ and form the nm t h matrix element 
of the partially coherent density operator as 

1 


< n,p(z;c)m > = e ^ ^ 1+< ^yT 


+ c(l + l/c)( m+n / 2 ) 

l/ 2 / _ \ ri-m / U |2 


S) fe) *-(<&) • <*> 
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where P, m (r) is the polynomial given by 

m I 

P m ( r \ = r k 

1 V ’ z^k\{l + k)\{m-k)\ 

The thermal or noise density can be found from the c — » 0 limit of the previous equation 
and is given by the familiar expression 


< n, p(z ; 0)m >= e 




(29) 


VVdra! 

By calculating another Gaussian integral, an overlap of two partially coherent states is 


found to be 


Tr[p(z 1 -,c 1 )p(z 2 ]c 2 )] = 
Clearly as c\ — » c 2 = c, eq. (30) reduces to 


exp 


1+Ci+Ca 


(1 + C\ + c 2 ) 


(30) 


Tr[p(0;c)p(0;c)] = y 


■I - 2c 


(31) 


Remark: The non-zero part of the entropy 5(c) = 5(z; c) and the non-unit part of 
Tr ( p{z , c)) 2 give a measure of the nagnitude of the departure from a pure coherent state. 
The generalization to M modes, with M a positive, finite integer is straightforward. Let 
(a,z) be M x 1 matrices, let (a*, z*) be 1 x M matrices, 1 the M x M unit matrix where 
k = 1,2, ...,M labels which mode. Let B be a given non-singular, Hermitian, M x M , 
positive, covariance matrix and express the M-mode Fock state as 

M (a*) n * 


|n!,...,n m >= TT — =-|0 > , 
t=l vW 


(32) 


where is the occupation number of the k th - mode. The M-mode thermal density operator 
is given by 


m(B) = <fcf(l + B) (33) 

-imSW z ' e ' x '' e '^ x><A ’ (34) 
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where x = (xi, , , . , x m) is the Af-component, complex, coherent state amplitude. The 
entropy of an A/-mode phonon packet is given by 

S(B) = k B Tr[(l + B)ln(l + B)-BlnB]=S(x\B) . (35) 

for both the incoherent and the partially coherent cases. The pair correlation function g 1 ' 2 ' 1 

for an M-mode partially coherent system is 

M) (z . B) = Tr[B_B + 2z^_Bz] 

9 K ' } (Tr(B) + |z| 2 ) 2 

which has the limits 

9 « 2 >(0,B) = 1 (37a) 

and 

9 < 2 >(z,0) = 0 (376) 



To find the partition function for the partially coherent system one needs the mode energy 
e* for the A;* ft -mode which is real 


5 


(38) 


and the A<f-mode vector 


e = 


/ex X 

\e_M ) 

which is a vector-valued, intensive variable. The mode energies for the k th mode is 


(39) 


E k = n k e k - El 


(40) 


where n k is the occupation number for the k th mode and 


£= (E u ...,E m ) 


(41) 


is a real, vector-valued, extensive variable dual to eq. (39). In this more general case the 
density operator p{z\ B) which maximizes the entropy 


S - —kgTrlp Inp ] 
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is subject to the constraints 


(42») 
(426) 
(42c) 

( 42rf) 

Let (3 be an m-vector with the value j3 = 1/kgT for each component and express the 
given covariance matrix for bosons as 

B = B(f3) = (e~ 0 e - l)” 1 = {e~ fiE - l)" 1 . (43) 

( A similar argument for Fermi-Dirac particles would replace —1 by +1 in eq. (43) but 
this is not needed here.) The partition function Z(/3,V ) is now 

Z(l3,V) = d e tf _ 1 e _ IIE 'j (44) 


and 


Tr(p) = | , 

Tr(pa) = z , 

Tr(pa*) = z * , 

Tr(pa* e a) = z* e z+ e 


at thermal equilibrium and V is the volume of the fluid plus radiation. From eq. (44) all 
of the equilibrium thermodynamic quantities can be calculated, for example the average 
energy S is 


££Ve-^ r 

£ = — 

z 


dln(Z) 

d(3 


(45 a) 


when classically 

Z(/3,V) = J2^ Er (45 b) 

r 


and in quantum statistical mechanics a Trace over matrix elements of e is taken. The 
pressure p is 


1 dln(Z) 
£ dV 


(45c) 
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and the (reversible) differential work dW is 


dW = —pdV 


(45 d) 


and the previous equation can be used to obtain 

13 dV 

The entropy of eq. (10) can also be written as 


dV 


S = k B [ln{Z) + (d£\ 


(45e) 


and the Helmholtz free energy is 

fXAV) = -ifci[Z(/J,V)] . (45/) 

The coherent state density operator can be given as 


p(z) = det 


(1 _ e -pE^ e -(c*-* m )0E(a-z) 


(46) 


Near equilibrium, a Kubo linear response theory 21 can be established by studying small 
deviations from equilibrium where the fluctuations will be equal to the dissipations. This 
exercise will be left to a future project. 


3. Conclusions and Outlook 

A two-component thermodynamics was formulated for a fluid in thermal equilibrium 
with a reservoir radiated by a coherent state of phonons. One future project will be to 
derive the Kubo fluctuation-dissipation theorem for this system, another will be to compare 
and contrast these results with both the Langevin equation and the stochastic quantization 
approaches. These future studies should illuminate (or simplify) the lattice Monte Carlo 
methods. 

Much remains to learned from harmonic, SHO , systems. 
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